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ABSTRACT 


An analytical model of a jet injected normally from a flat 
plate into a uniform crossing flow was formulated to provide 
a simplified method of predicting the interference effects 
arising from the comolex flow fields induced by ship bow 
thrusters. This model was an extension of previous work 
based upon a description of the jet as a series of distributed 
vortices. The analysis takes into account the position of 
the effective source of the jet and the blockage due to the 
presence of the jet in the crossflow. For representative 
jet-to-crossflow velocity ratios, the flow field and pressure 
distributions were calculated utilizing different combinations 
of effective source position and blockage. The accuracy of 
the model was evaluated by comparison with the available 
experimental data. Although good agreement was achieved for 
large portions of the interaction field, several regions were 


identified requiring further analytical description. 
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fet NT RODUCT ION 


In recent years there has been an increase in the instal- 
lation and use of hull mounted, tunneled bow thrusters to 
improve ship maneuverability in restricted situations, such 
as dynamic position holding and confined waters docking. 

This increase has created renewed interest in the development 
of an ability to predict bow thruster performance. To date, 
there has been relatively little pre-installation design 
consideration regarding precise power requirements or ultimate 
operational performance. The limited testing that has been 
conducted / 1-4/* has resulted in a recurring major problem, 
namely, radical variations in effective turning moment with 
minor changes in ship speed. Chislett and Bjorheden / 4 7 
noted that "an area of low pressure is created downstream 

of the discharging jet producing a resultant suction force 
that has a shifting center of action with changes in shiv 
speed," thereby altering the effective moment produced by 

the thruster. Tne ability to predict the magnitude of this 
resultant suction force and its center of action would permit 
designers to accurately determine powering and control re- 
quirements for thrusters. 

This problem of radical variations in effective side 
forces and turning moments with only minor changes in ship 


speed may be better understood by considering the general 


Numbers in brackets refer to references listed. 


equations of motion of a ship moving in the horizontal 


plane of the sea. From / 5 / these equations are: 


m(u - rv) =X 
m(v - ru) = Y 
Ir=N 


Z 
where m= mass of the ship 

u = velocity in the x-direction 

v = velocity in the y-direction 

r = angular velocity 

I = mass moment of inertia 

X = excitation force(s) in x-direction 

Y = excitation force(s) in y-direction 

N = excitation torcue(s) 

(dots indicate differentiation with resnvect to time) 
To minimize the complexity of the following brief analysis, 
all forces and torques will be referred to the shio's center 
of gravity. See Fig. l. 

When these equations are transformed into their linearized 
version / 5 /, they become 
= <u Re? tk) ul = Xx 


<= V = © . = = = ee = 
YY + (m Yelv (Y mu.)xr Yer 4 
= Nov - Nev = Nr +. (I, - Ny)r = N 
where the subscripts u, u, Jae v, 12 5 r, denote differentiation 
with respect to the variable subscripts, and uy is constant 


reference surge velocity. The left sides of these equations 


explicitly represent the velocity and acceleration dependent 
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reactions due to drag, added mass, added inertia and other 
coupled motions; while the right sides implicitly represent 
excitation force(s) and moment(s), which are in this analysis 
due to the thruster. Since this is a linearized formulation 
utilizing small angle perturbations, the excitation force in 
the x-direction 1s negligible compared to that in the 
y-direction: 
mise 
therefore X=0. The y-direction excitation force can be 
expressed as 
Yea Yee by 


7 a 
the ideal jet thrust produced by the 


(1) 


here Y. = m.U. 
oe ares 


thruster, and Y; 1s the induced force created by the inter- 
action of the thruster discharge with the flow past the ship. 


Utilizing experimental data from / 4 / in the range, 


Orel < 5 SUmeo, Y; can be approximated by a linear function 
J 


of ship speed 


3 


g— = 0.22 + k(t) (2) 
j U5 


where -3.0 < k < -1.5. Substituting a nominal value of 


k = -2.0 and Y. = mU, into equation (2) gives 


a> U 
Y; = 0, [0.22 ~ 2.0 | (3) 


Substituting equation (3) into equation (1) results in 


> U 
Y= BU. nee - 2.0(— 
HO J ‘3, 
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which explicitly couples the y-direction motion to the 
x-direction motion, normally uncoupled in the linear 
approximation. In addition to the coupling, in the relevant 
velocity ratio range, 0.1 < 7 < 0.25, there 1S an appreciable 
reduction in the net thruster force due to motion in the 
x-direction. Furthermore, the moment produced by the thruster 
is a function of the net effective thrust, Y, and a character- 
istic moment arm, x, ? 

Ny = x, ¥ 
As reported by Refs. /1-3/, the degradation of net effective 
thrust is accompanied by a shift in the center of action of 
the induced force, Yes which results in a change in the char- 
acteristic moment arm. This shift in the thruster's net 
effective thrust center of action is a direct consequence of 
the coupling of x-direction motion with turning moment. 

This simplified analysis, which neglected all excitations 
other phan the thruster, points out the complex, coupled mo- 
tions arising from the interaction of the thruster discharge 
and the flow past the ship. To fully understand this highly 
non-linear, three-dimensional flow situation extensive 
analytical and experimental work is required. This thesis 
will attempt to provide some understanding of the interaction 
of a jet discharging into a uniform cross-flow in order to 


develop some ability to predict the degradation of net effective 


EnNrUSt. 


v2 





II. BACKGROUND 


A survey of the literature pertaining to the installation 
and operation of bow thrusters revealed a surprising lack of 
experimental data, considering the number of units currently 
installed. English eae Stuntz and Taylor a and 
Chislett and Bjorheden / 4/7 present some data from model tests 
conducted at various facilities, but in these works there is 
no formalized recording of full scale test or actual instal- 
lation test data. However, further investigation led to a 
related field: Vertical or Short Take-Off and Landing Aircraft 
(V/STOL) research. The concept of using air eae positioned 
in aircraft wings to provide vertical take-off and landing 
Capability correlates extremely well with flow patterns pro- 
duced by a bow thruster. As reported by Bradbury and Wood 
/ 6_/ for incompressible flow, the jet (thruster discharge) 
path and induced flow are dependent mainly on the momentum 
mix ratio, 950, °/pU", and independent of Reynolds number. 
Margason Va 7, and Gordier ae, also reported that the 
effective velocity ratio is the predominate characteristic 
in determining the path of the jet. Therefore, until more 
extensive ship or model data have been collected, the re- 
sults of V/STOL research are used as a basis for predictive 
theories. 

Although there has been a good deal of research conducted 
in the area of air jets in cross-flows, the majority of the 


worx has been concerned with analytically or empirically 
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defining the jet trajectory and jet cross-section geometry. 
The interference effects of the interacting flows have re- 
ceived relatively little attention. A fluid jet injected 
into a crossing stream has previously been described as a 
turbulent, three-dimensional, highly non-linear flow problem 
which, even with considerable simplification, requires ex- 
tensive computer time to numerically solve the appropriate 
form of the Navier-Stokes equations / 9 7. Therefore, early 
investigators have sought simplified methods for determining 
the interaction effects, such as pressure and velocity dis- 
tributions, and have realized the necessity of knowing the 
geometry of the jet: trajectory and shape. 

Experimental data form the basis for most trajectory 
formulations. Either pressure or velocity measurements have 
been taken in the flow field encomvassing the jet-crossflow 
interaction region and curves fit through points of maximum 
pressure or velocity. These curves were then compared with 
flow visualizations. In some instances, semi-empirical for- 
mulations were derived from the conservation laws combined 
with experimentally determined constants. Abramovich / 107 
qualitatively described the turbulent jet in a deflecting 
flow and presented some empirical methods for predicting 
the trajectory of the deflected jet. Jordinson /11_/ pre- 
sented trajectory data by recording and plotting contours 
of total pressure coefficients. Keffer and Baines / 12 _/ 
presented experimental results from which it was determined 


that for various jet strengths, the jet trajectories could 
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be represented by a single function. Gordier / 3 7, the 
only reported source of water jet-water crossflow experiments, 
presented an empirical jet trajectory formulation based on 
a curve passed through experimental points of maximum total 
pressure. Sucec and Bowley / 9 7 formulated an analytical 
expression for the jet trajectory utilizing previous experi- 
mental information and the assumption that the distributed 
pressure force and entrained momentum flux could be approxi- 
mated by an aerodynamic drag force. Margason a utilized 
flow visualization techniques and pressure measurements to 
formulate a trajectory equation. In addition, his paper also 
reviewed and compared the results of a number of other 
studies concerned with trajectory prediction. All the above 
predictive formulations produced trajectories within a range 
of uncertainty that can be attributed to expected experimental 
error, differences in test procedures and facilities, and in 
the case of analytical formulations, simplifying assumptions. 
In view of this, it was concluded that any one of them re- 
flected the current state of the art in trajectory prediction. 
The ability to predict the jet trajectory is of definite 
importance, but the primary objective of this jet-crossflow 
research is to provide insight into the changes, in the other- 
wise uniform flow, created by injecting a jet. Experiments 
conducted by Bradbury and Wood / 6 7, Vogler / 13_/, Fearn 
and Weston / 14 7, McMahon and Mosher / 15_/7, Kamotani and 
Greber / 16_/ and others, provide measured values of the 


pressure distribution on the surface surrounding the jet 


BL) 





orifice which have been used to compare the accuracy of 
predictive models. During the past few years, a number of 
approximate methods to voredict the pressure distribution has 
been developed: Wooler /17,18/7, Wu and Wright [SV Adler 

and Baron aeOP and Schmitt ON Again, because of the com- 
plex nature of the flow, all these models have utilized some 
degree of empiricism. Either the model is formulated using 
experimentally obtained trajectory expressions or the model 

is formulated using analytical expressions whose coefficients 
are selected to provide best fit with experimental data. 

The majority of the models are based on integral techniques 
with simplifying assumptions, such as (i) the representation 
of the entire flow as two-dimensional Le (11) the external 
flow is irrotaticnal, uncompressible and steady-state ERZO7 ei 
(111) the flow is inviscid except that viscosity is the 
mechanism that leads to entrainment /18 7. 

With the aid of these assumptions, the following approaches 
were made. Adler and Baron formulated their model by integrat- 
ing momentum equations to describe the jet mixing field, without 
using empirical trajectory data. However, numerous parameters 
were derived using empirical correlations from prior research. 
Their results provide satisfactory agreement with experiment, 
but from the quantity of empiricism involved, one might ex- 
pect such agreement. Wooler all 7 followed by Wu and Wright 
VS 7 utilized blockage-sink representations where entrain- 
ment of crossflow fluid was handled analytically with suitably 


chosen coefficients. Again the results had satisfactory 
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agreement with experiment. However, an earlier effort by 
Wooler / 18_/7 minimized empiricism to the extent that only 
an experimentally determined jet trajectory was necessary 
to complete the model. By thus restricting the number of 
empirical parameters, this model has considerable anpeal 
from an engineering point of view. Since visualization of 
actual jet-crossflow situations has shown that the jet is 
deflected in the crossflow direction and forms two contra- 
rotating trailing vortices, this interference model was 
based on the representation of the jet by a distribution of 
vorticity. Incorporating arguments from the aerodynamic 
theory of lifting bodies this vorticity distribution was 
quantitatively described along the experimentally determined 
trajectory. By virtue of minimal empiricism and therefore 
greater appeal, this model was selected as the basis for 
continued investigation in developing a better predictive 


E@O |. 
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IIIT. REVIEW OF WOOLER'S VORTEX MODEL 


Before any attempt to develop an improved model was 
made, a detailed re-derivation of the basis formulation was 
performed in order to: (1) gain a better understanding of 
the logic of formulation and (2) confirm the reported results. 
A review of the distributed vortex method follows. 

One area of interest was the geometry of the problem. 


The empirical trajectory equation used was 
x/d = B [ cosh (2/Bd) - 1 | (4) 


where x iS measured in the direction of the mainstream and 

Z 1S measured in the direction of the exiting jet. The 
coefficient B was empirically determined from Jordinson /11/ 
and is equal to 0.19 (U,/u)*, where o is the jet velocity 
and U ie the mainstream velocity. To further describe the 
problemi{(see Fig. 2), a system of natural coordinates 
attached to the jet is adooted. Relative to the fixed 
coordinate system and an arbitrary point Xp r¥p| on the 


surface surrounding the jet orifice, the natural coordinates 


are given by 


gE= - | zsina af (x-x,,) cosa | 
n = a (5) 
C= 


ZCOSa - (X-X,,) sina 


where € is the coordinate tangent to the jet; ¢, the coordi- 


nate normal to the jet in the direction of the center of 
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Gurvature; n, the coordinate perpendicular to € and ¢f ; 
and a , the angle between the x and € directions. 

Another extremely important part of Wooler's formulation 
was the calculation of the distribution of vorticity within 
the jet. This distribution of vorticity was determined by 
adopting methods commonly found in the aerodynamic theory of 
lifting surfaces / 227. According to this theory, regions of 
flow external to the jet are taken to be irrotational and the 
deflection of the jet is due to a purely inviscid mechanism. 
This inviscid mechanism 1s expressed as a balance between the 
forces due to pressure differences across an element of the 
jet and the centrifugal forces associated with the jet curva- 
ture. These pressure differences define corresponding velocity 
differences according to Bernoulli's Theorem. These resulting 
velocity differences, in turn, are related to circulation 
according to Kelvin's definition of circulation. Therefore, 
with the above assumptions, the circulation is related to 


jet curvature, and this relationship is 


Z 
u 
rT soared 7 
ou - 7-4(74) Ho 


where [T is the circulation per unit length, ds, along the 
jet; R, the local radius of curvature of the jet; and d, the 
diameter of the jet orifice. Utilizing the empirical tra- 
jectory expression / Eq. (4)/ in the arc length derivative 


results in 


ae ares az (7) 
R Bd cosh (z/Ba) 


Bs, 





imeorporating Eq. (7) into Eq. (6) and integrating over an 
element of jet extending from Z4 to Z5 results in the follow- 


ing expression for the circulation of a jet element 


2 
a Z Z 
va (73) [ tan” 2 2 73a - tan € ‘ea | 


“A 
il 
bof a 


: (8) 
In order to determine the flow field disturbance on the plate, 
the jet was divided into a number of vortex elements whose 
strengths were calculated using Eq. (8). These elements are 
the origins of characteristic horseshoe vortex systems where 
the cross member is bound in the jet and the trailing sides 
are, according to Wooler's assumption, tangent to the jet and 
separated by a distance, d, the diameter of the jet orifice 
(Fig. 3). Errors associated with taking the trailing vortices 
tangent to the jet instead of containing them within the jet 
as in the actual flow situation were assumed by Wooler to have 
negligible effect on the velocity field along the plate. Each 
of these horseshoe vortex systems produces an interference 
velocity on the surface surrounding the jet orifice. To de- 
termine this interference velocity at any arbitrary point on 


the surface, the Law of Biot and Savart / 237 


> => 
- Ee K_ ds x xr 
q LZ = 


is applied to each bound vortex and its associated trailing 
vortices. For the bound vortex, - di2 <L<§ oe and for the 
trailing vortices, zs L<o , were the intesration limits. 
These integrations produced closed-form solutions for the 


interference velocities due to individual jet elements. The 
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total interference velocity components due to the entire jet 
are determined as the sum of elemental contributions and 


are given by 


/Bd z,/Bd 


u a Oe -l/e i+] - e 
o = ud » tan ee x[(-w, tw)-w,) sina, tu, cosa, ] 
CZ. + z.)/Bd 
: l+e* "itl 1 
1=l1 
, . U3 _ 441/34 _ Zs ,/3d 
ri Gp ate so ee 
i=l - 


where N is the number of jet elements and Zs and Z are the 


i+l 
endpoints of the general element. The parameters u, v and 
w are defined as geometric coefficients resulting from the 
integration along the jet. The subscripts l, 2, 3 refer to 


the bound, left and right (as viewed from upstream) vortex 


contributions respectively. They are given by 


= Se 
i 
Pack gon? +074 (pte)? Vf *4e 7+ in- 4)? 





C g 
vy, =e SOURS 
f p*t(n= 8)? , [e2+¢7+ (n- 5)? 
G g 


v, = Qo 
5 tent 8)? ete (ne? 
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nb n+ 1 - + 
_- : 
. SS 
E +e \e “i rc? oa 1) ¢ g? ni r? ener x.) 4 
wy = O04 ie a 
co + ( ne)? eo to + in - 4)? 


rf 


Ww, = Nn + 3 ja & 

2 2 2 

c7 + (mn + 5) ee tees) 
and 6, & , and n are given by Eq. (5) for each point on the 
surface. After calculation of the interference velocity 
vector at a sufficient number of voints on the surface, the 


pressure coefficinet at each point, Co? is given by 
P- PL 
2p are 
= =jl- 
p pe q /U 
2 oU 
—— Bs au +(— ae The double u and v 


components are the result of the use of an image system to 


ec 


where ue 


establish the surface as a solid boundary thereby negating 
the w component of velocity in this plane. Upon substitution 


and simplification, one has 


[eG + | 8 


After completion of the review and the incorporation 


Q 
ll 


of several necessary corrections to the published work /18 /, 
the entire formulation was coded for computer solution. The 


resulting plots of pressure contours for representative 
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velocity ratios (Figs. 4,5,6) show fair agreement with the 
experiments of Bradbury and Wood / 6 _/. However, these 
plots also point out serious deficiencies in the model. In 
the far field, in the arc 35° < @ < 145° (68 being measured 
counter clockwise from the ray extending downstream from 
the jet origin), the method represents the jet interference 
effects well, but in both the upstream and downstream 
portions of the surface, there is a definite lack of agree- 
ment. In the upstream area, the blockage effect due to the 
presence of the jet in the uniform flow is not predicted, 
while in the downstream regions, additional wake effects 
have not been taken into account. However, in view of the 
large area of a agreement with experiment, Wooler's 
distributed vorticity model was selected as the foundation 
for further investigation. It is apparent that this model 
gives a fair representation of the actual contra-rotating 


vortices associated with the jet-crossflow interaction. 
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IV. EXTENSION OF WOOLER MODEL 


Upon successful confirmation and correction of Wooler's 
distributed vortex model and as a step towards developing an 
improved model, it was felt necessary to verify the trailing 
vortex assumption, that these vortices are tangential to the 
jet rather than contained in the jet. Retaining the basic 
theoretical concevts for calculating the incremental circula- 
tion of a finite length of the jet, containment of the trailing 
vortices was approximated by successive conforming discretized 
vortex elements. These vortex elements were formed by restrict- 
ing the length of the trailing vortices to the linear distance 
between Zy and Zor the arbitrarily-chosen end points of the 
finite length of jet, ds. In order to satisfy the Helmholtz 
vortex theorem, that a vortex cannot end in tne fluid, the 
above horseshoe vortex system was closed by connecting the 
trailing arms with another bound, contra-rotating vortex. 

(See Fig. 7) By making this alteration, vortex loops of in- 
creasing incremental strength, as given by Eq. (8), can be 
made to conform to the curvature of the jet, as in the actual 
flow situation. This resulting vortex pattern is analogous 
to that created by a continuously accelerating wing / 23 7. 
As a wing accelerates, the strength of the produced vortices 
correspondingly increase, a situation very similar to the 
change in elemental vortex strength associated with the jet 


curvature. 
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Using this contained-vortex system, the interference 
velocity due to the jet at any arbitrary voint on the plate 
was determined by applying the Biot-Savart Law (See 
Appendix A). In this application, it was necessary to 
approximate the trailing vortices as straight-line segments. 
However, as the sizes of tne jet elements become small, the 
linear approximation can be made to conform to the curva- 
ture with an exactness that is only limited by numerical prac- 
ticalities. The total interference components due to the 
entire jet result from the summation of elemental contribu- 


tions and are given by 


u 1 U5 ‘= : ; : 
7T7~ a\q » mw sina, te,sinay, + (-wytwg) sina, tu ycosaymuzcosa, se 
i=l ) 
1 Z Z 
a m Kel bs k 
oi (zy he 
jee k+1 k 
k=1 = 
2N a Zz Z 
vi l( 4 vera vt x tan™? a 
G 8 \U 3 4 (Zr 44 4 Zz) 
a k=1 ee 


where N is the number of jet elements and Z1 and Zy4, are 
the endpoints of the general element. Again, the parameters 
u, v, and w are defined as geometric coefficients resulting 
from applying the Biot-Savart Law around the element. The 

subscripts 1,2,3 and 4 refer to the forward and after bound 


and left and right trailing vortex contributions, respectively. 


These are given by 
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o n+ now 4 
u = ee a 
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eu oe rma jen Pea Gh - hye fi 
ca Ct (n = 35) 532 +o, + (n - 4) eBl +o, «+ (n = 35) 
SB SB? Bl 
a> oD uaa 2 
can Ct (n + &)? E32 +o, + (n + 35) eBl +o, + (n + +5) 
n+ 3 no- 4 
a joe eee 
cues t (n + 35) by + 5, + (n = 4) 
50 n+ 3 n- 
Ww, = = 
gt cen? + cy + in te? Ven? + 052 + n= 8D? 
=a bl Bl 
w., = = 
3 re _— 2 2 2 2 
by = 3) vez, + 2) Si ee) 


n + ¥ >B2 "Bl 
es ae St Piney? Vet, + cae - in + 4D” 
where the various &'s, t's and n's are given by Equations 
(Al, A2, A3) of Appendix A. With these interference velocity 


components at numerous points on the plate, the pressure 


coefficient, Eq. (9), was determined for selected points 


26 





on the plate. The resulting contours of constant pressure 
coefficient are shown in Figs. 8, 9, and 10. 

These resulting plots show a somewhat different agree- 
ment with the experimental data / 6 7 than did the results 
of Wooler's model; however, the discrepancy can be directly 
attributed to the manner in which the trailing vortices are 
treated. In Wooler's model, the total interference velocity 
due to the trailing vortices results from the cumulative el- 
emental vortices acting at different distances relative to 
a point on the surface, while in the conforming vortex model 
the trailing vortex contribution results from the cumulative 
elemental vortex strength acting at a single radius, (See 
Fig. 11). Therefore, Wooler's assumption of treating the 


trailing vortices as being tangent to the trajectory intro- 


duced an error that under some conditions tended to improve 
agreement, where the conforming vortex model removed this 
error at the expense of some loss in agreement, at least at 
the larger values of U;/U. However, there have been 
studies VE 22, that present the concept of an effective 
origin of the jet vortex system - a region where the turbu- 
lent mixing processes extend across the entire jet causing 
definite deflection of the jet and establishment of the contra- 
rotating vortices. Schmidt / 21/7 formulated an empirical 
expression for the position of the effective source as a 
function of the jet-to-mainstream velocity ratio given by 


Z = (15 ava) "> 
e 


where oO = Wee - Applying this effective source expression 
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to the conforming vortex model produces better agreement 
with experiment (Fig. 12). 

Comparison of the conforming vortex model, including an 
effective source, with exverimental data shows the same, 
if not better, far-field agreement with an improvement in 
the predicted pressure distribution in upstream and downstream 
regions. Therefore, this model which more closely represents 
the actual flow situation, replaced Wooler's model as the 


basis for further investigation and improvement, 
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V. IMPROVEMENT OF CONFORMING VORTEX MODEL 


Evaluation of the conforming vortex model results in- 
dicated the continued lack of agreement with experiment in 
two distinct areas: uostream and wake. In order to improve 
the agreement in these areas, it was felt that additional 
corrections could be formulated from a qualitative analysis 
of the flow phenomenon. The discrepancies in the upstream 
area can be, at least partially, attributed to the blockage 
of the mainstream flow by the presence of the jet, while the 
discrepancies in the wake area are caused by turbulent flow 
separation, vortex shedding and turbulent entrainment. In 
view of the present uncertainty concerning the exact character 
of the wake area, it was decided that the development of 
corrections for the wake must await future investigation and 
experimentation. Therefore, the upstream blockage correction 
was undertaken. 

Near the point of injection, the physical appearance of 
a circular jet injected into a crossing flow is that of a 
cylinder. Interpreting the jet as a cylinder in an irrotational 
Flow, potential flow theory predicts perturbations to the 


uniform flow that are given by 


Z 
= == cos 28 
(10) 
Z 
Vv 
U =-55 Sin 28 
r 
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where r and 6 are cylindrical coordinates of points in the 
flow field and ais the radius of the cylinder. These 
velocity perturbations caused by blockage can be combined 
with the interference velocity components derived from the 
vorticity formulation to adjust the coefficient of pressure 
(Eq. 9) at arbitrary voints on the surface. The initial 

(and most simple) assumotion was to make the diameter of the 
blockage cylinder the same as the diameter of the jet orifice. 
The resulting pressure distribution indicated an excessive 
blockage that is not present in the real flow. The actual 
flow does not produce a solid interface, as would a cylinder, 
but instead an entraining interface surrounding a core of 
uniform velocity approximately equal to Vis the jet exit 
velocity. To better describe this, a potential cylinder of 
diameter less than the diameter of the jet was introduced. 
Using many combinations of cylinder diameters with different 
Jet -to-mainstream velocity ratios, an expression for a in 

Eq. (10) was formulated in terms of the velocity ratio. 


a= 0.96 vo 


Applying this formulation to the blockage-cylinder perturba- 
tion velocity calculations results in pressure contour plots 
(Fig. 13,14,15) that show improved agreement with experiment 
in the far field and upstream areas for jet-to-mainstream 
velocity ratios, 4 < < 8. However, as the velocity ratio 
increases, the model results and experiment diverge. The 
apparent cause of this disagreement can be resolved by 


qualitatively examining the actual flow. With increasing 
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velocity ratios, the entrainment of crossflow fluid by the 
jet becomes the more dominant interference factor, while the 
blockage factor is relatively less important. Wu and Wright 
/19 7 have concluded that "the amount of crosswind fluid en- 
trained increases with increasing jet speed" when the cross- 
wind speed is held constant. In addition, this increased 
entrainment appears to be especially important in the down- 
Steam wake region and causes a further reduction in the 
blockage effect as compared to that due to a solid cylinder 
in a crossflow. The resulting effect of the increased rate 
of entrainment at higher jet speeds is larger interference 
velocity perturbations in the flow field surrounding the jet 
orifice near the plate. In view of this, the lack of agree- 
ment at higher velocity ratios can be explained while also 
pointing out the need for additional correction(s) to account 
for the change in the relative importance of blockage and 


entrainment with changes in velocity ratio. 
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VI. CURRENT MODEL CAPABILITY 


Realizing the limitations of the current analytical 
model, it is nevertheless appropriate to investigate its 
predictive capapility. Figure 16 shows the results of the 
analysis expressed as the force induced on a flat plate due 
to a jet injected at right angles to a crossflow. The force 
1s expressed as a fraction of the ideal jet thrust and the 
negative values indicate suction forces (in a direction on=~ 
posite to the ideal jet thrust). Correlation with experiment 
/ 4 7 shows a difference in the value of k in Eq. (2). The 
value of k = ~4.4 resulting from this analysis, although not 
in the predicted range, ~3.0 $k S$ -1.5, still displays the 
correct trend and provides encouragement as to the potential 
success of a model based upon the present methods. The main 
difference in the approximation constant results from the 
lack of proper accounting for the entire entrainment effect, 
a factor which progressively degrades the model at higher 
velocity ratios. However, as an initial design approximation 
for low velocity bow thrusters active at relatively high 
ship speeds, the conforming vortex model presented here pro- 
vides a satisfactory prediction of the interference effects 
-- a step toward improving the accuracy of vowering and con- 


trol determinations. 


32 


VII. CONCLUSIONS 


The conforming vortex model, modified with corrections 
based on experimental observations, is presented as an 
initial formulation of a simple yet physically consistent 
representation of jet-crossflow interference effects. Pres- 
sure distributions in the far-field are adequately predicted 
and the calculated results obtained with the model provide 
reasonably good agreement with experiment for the lower 
jet-to-crossflow velocity ratios. Agreement is less than 
satisfactory for the higher velocity ratios. However, the 
model does indicate the relative importance of different 
flow phenomena with changing velocity ratios. The most evi- 
dent weakness of the model is the lack of an adequate method 
for including entrainment and blockage effects near to the 
origin of the jet at higher velocity ratios. By restricting 
the range of application, the conforming vortex model demon- 
strates sufficient agreement to be utilized as a first 


approximation technique for designing bow thrusters. 
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VIII. RECOMMENDATIONS 


Further experiments should be conducted to provide more 
correlation data, especially in the case of water jet into 
water crossflow. Empirical constants utilized in the current 
model and those that will become necessary when trying to 
model increased entrainment can be refined when a larger data 
base is available. The encouraging results of this method 
of analyzing the flow field should not be overlooked, and 
additional corrections, as indicated in this analysis, should 


be made to improve its accuracy. 
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Figure 1. Sketch of bow thruster-crossflow interaction. 
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Figure 4. Comparison Of Pressure Coefficient Contours 
Experiment ---/ 6 _/, Wooler 


Velocity Ratio = 4 
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Figure 5. Comparison Of Pressure Coefficient Contours 
Experiment ---/ 6_/, Wooler 
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Fiqure 6. Comparison Of Pressure Coefficient Contours 
Experiment ---/ 6 7, Wooler 
oe = 11.3 
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Figure 7. Conforming Vortex Loop System 
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Figure 8, Comparison Of Pueege rem coeiticzent Contours 
Experiment ---/ 6 7, Conforming Vortex Model 
Mas = 4.0. 
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Figure 9. Comparison Of Pressure Coefficient Contours 
Experiment ---/ 6_/7, Conforming Vortex Model 
U./U = 8.0 
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Figure 10. Comparison Of Pressure Coefficient Contours 
Experiment ---/ 6 /, Conforming Vortex Model 
yu = 11.3 
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Figure a. Wooler Model 





Conforming Vortex Model 


Figure b. 


Comparison Of Wooler (a) And Conforming Vortex 
Model (b) Effective Radii For Induced Velocities 


Figure ll. 
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Figure 12. Comparison Of Pressure Coefficient Contours 
Experiment ---/ 6 7, Conforming Vortex Model 
With Effective Source Correction—-U./U = 8.0 
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Comparison Of Pressure Coefficient Contours 
Experiment ---/ 6 /, Fully Corrected Conforming 
Vortex Model——.U, /U = 4.0 
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Figure 14, Comparison Of Pressure Coefficient Contours 
Experiment ---/ 6 7, Fully Corrected Conforming 
Vortex Model—.U,/U 76.0 
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Figure 15. Comparison Of Pressure Coefficient Contours 
Experiment ---/ 6 /, Fully Corrected Conforming 
Vortex Model—,U,/U = TE 
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Variation In The Induced Force Due To Jet-Crossflow 
Interaction 
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APPENDIX A 


FORMULATION OF INTERFERENCE VELOCITY COMPONENTS 


The formulation of the interference velocity due to the 
vortex loov system of the conforming vortex model is given 
below. The formulation requires the use of three separate 
coordinate systems that are functions of the endpoints of a 
jet element and the position of an arbitrary point on the 


plate. From Fig. (Al), these are given by: 


Si ei [zysina, a (x) — x) cosa. | (Al) 
= Yp 
o) = ZjCOSsa, - (x4 - x) sina, 


for the upstream bound vortex 


E> = «= | z,sina, + (X5 - x) cosa.| (A2) 
oem tn 
So = 2,c08a, - (x, - x) sina. 


for the downstream bound vortex, and 


SBi = - |z,sina, + (x; - x) cosa,,| 
n = (A3) 
as 
Ca OF «O2yCcosa, - (x; = X,) sina, 


ot 








for the trailing vortices, where aS and Y,, are the 


coordinates of an arbitrary point on the plate and x Z 


ee 
X5,Z5 are the coordinates of the endpoints of an element 


1 


of jet. The angles Ayr Go and a, are the angles between the 


B 


x-direction and respective &'s and are given by 


a, = tan’ + |esch (z,/B)} 
a, = tan™> [esch (2.,/B) 
a, = tan? ze - Z,)/ (x, - x1) 


To determine the induced velocity at an arbitrary point, 


the law of Biot and Savart 


sing 2 
= / 5 as (24) 
r 
L 


must be individually applied to both the bound and trailing 


iF 





lq| = 


a 


vorticies. Due to similarities in the integrations, only 
one bound and one trailing vortex formulation will be pre- 
sented and necessary changes applied for the other vortices. 
The circulation of a finite amount of jet, given by / Eq. 


( 8)7 
2 
U. ‘ z./Bd 7 z,./Bd 
et va( 21) lean? 9”? Means et 
is the incremental contribution which is added to each 


successive element, such that the total circulation for 


a given element is given by 
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n 


a 
U. 7 Z. /Bd S z, /Ba 
r - $a (qi) | Jean7h e# Oa ce a 
n U 


(AS) 


Therefore, the induced velocity due to any of the vortex 
filaments associated with the n-th element, 1 sn<sQN, 
where N is the number of jet elements, can be determined 


from 


lq oe { sine ds (A6) 
‘ r 

where L is the length of the filament, 8 is the angle 

between the vortex filament and the radius, r, to an 

footehary POlnt on the plate. 


For the forward bound vortex (See Fig. (A2)) 


Se 2 


sing = —_—_,—— 
VE,° ae on ete s)* 
and 
i 
me = Ei + a + (n + s)* 


Substituting these equations into Eq. (A6) gives 





laa | 2 Vey" + 84 d (A7) 
Baucus | jeu 2. 7 ale Tt : : 
nl at le, A cy eons) 3/2 


L 


the magnitude of the induced velocity at an arbitrary point 


on the plate due to the bound vortex of the n-th element. 


The x-direction component of this velocity is given by 
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eo = lai! lcosy cosa, - siny,sina,| (A8) 
I at 
where  —————— and siny, ee 


2 2 2 


Inserting Eq. (A7) into Eq. (A8) and applying the integration 


limits gives 






5 
r oy 
n 
—) = ea | ee, a ds cosa, — 
=k eer sme sl 
(A9) 
5 
/ BS 
a ds sina, 
as fe? + c,2 + m+ 80] 
Integration of Eq. (A9) results in 
es 
We = qn 3(—t ee ‘cosa, 
ae i - 1/61 eel ch )¢ lev fen 
= <a Daas = ee tine) 
51 Filet + cf + (nt)? fet + cf + (n-¥) 


(A10) 


the x-direction component of induced velocity due to the 
forward bound vortex. Similarly, the x-direction component 
of induced velocity due to the downstream bound vortex is 


given by 
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q.. = = a Pt Ss n- 3 \cosa, 
x2 T 
VE 5 ts aoe (n + 1)? VE oth, ene es 2) 
- ia e sina. 
Fe P45 24 (na 2 ne Zant ra 5 + n= Mas 
(All) 


The approximation of the bound vortices as being linear 
eliminates the y-direction component of induced velocity 
and since an image system is utilized to formulate the flat 
plate, it is unnecessary to calculate any z-direction 
component of induced velocity. - 

For the left trailing vortex filament, as viewed from 


upstream (See Fig. (A3)). 


ay eee? 2 
cn Ct (n - 3s) 


——— ————— 
ty 2 Z 2 

on Ct (n- 4)” +s 

r2 eae ee 1) 2 52 


Inserting these expressions into Eq. (A6) gives 








ye ne rae 5)? 
= c—- S 
Gn3 eB | : : ; (A12) 
= cn tCt ( n - 3)" +s 


the magnitude of the induced velocity at an arbitrary 
point on the plate due to the left trailing vortex of the 
n-th element. Resolving this velocity into x and y com- 


ponents gives 
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ao = Gn3 cosy, sina, 
and j (Al13) 
Iy3 In3} 5*7Y_ 
Wem *: SB 
0 and siny, = 
2 2 2 Z 
Sa Ot O(n 4s) Ea + (nga) 


Inserting Eq. (Al2) into Eq. (A13) and applying integration 


limits gives 


n- , 
q = 8 Cahn as sing 
Bo ne bait (n=) 2482 ] 97 Z B 


: (A14) 
r B2 on 
g 8 = rq J ds 
y Le? 
tals + (n=) “+s ita 


Integration of Eq. (Al4) yields 





m ee °Bl 
Wy 3 ~ 9" Th i; im) aa. > ae - SS aa sina, 
G ste CI seems Fos > / - 


ee a0 °B . *B2 ; "Bl 
y 3 4m 2 2 
C + (n+) 2 - Z Zz z = 2 2 
B Ven t(n-4) “+655 Vey + (nok) “+85, 
(A1L5) 


the x and y comnonents of induced velocity due to the left 
trailing vortex. Similarly the x and y components of the 
induced velocity due to the right trailing vortex are given 


by 
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n + 4 SB 


E g 
tT. + (n + 5 
B eG fa vee + (itil ren, 


B 


; Men oB °B2 7 SB1 
ep FD Ne 2 (nts) EZ Ven? (nth) 242, 


(Al6) 
Therefore, the total x and y components of induced velocity 


due to an entire element of jet are 


Gp * Ax *F Ix2 F ag + Fx 
(A17) 


q = q 


sya 3 eel 


y y4 
Inserting Eqs. (Al0), (All), (A15), and (A1l6) into Eq. (Al14) 


yields 


a 


qr yay cosa, - Us cosa,-w,Sina,+w,sina + (-w,+w,)sina,| 


n 
4am 
r (A18) 


where 
a et ont ne 
1 
+ ne 2 eS B 
51 PL IVE, “+0, °+ (nts) © VE, “+0, “+ (ns) 
<= 2 tae E ifn 
eo ae i. 9/5 Se 2 
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Substituting the expression for te from Eq. (A5) into Eqs. 


(Al8) and non-dimensionalizing yields 
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Figure Al. Coordinate Systems For Conforming Vortex Model 
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Figure A2. 





Geometry Of Bound Vortex Integration 
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Figure A3. Geometry Of Trailing Vortex Integration 
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C XXX HHH AANA HHI HH HAA KZ 
C X THIS PROGRAM SIMULATES A X 
c X FLUID JET INJECTEGS INTO X 
C X A UNITFECEM CROSSFELGh AND X 
C X FURTEER DETERMINES THE xX 
C X PRESSURE OISTRIALTICN ON X 
C X THE SURFACE FROM kKHICH xX 
‘ X THE JET 1S EJECTEC x 
¢ XXX XX AMX AAXX XXX ANAK XK AK K KK X 
C DEFINITICN CF TERMS 
c NIYFE — MODEL TYPE:C-CCNFORM 
C 1-CCNFORM SLOCKAGE 
C 2-RAISED CONFOR™ 
C 3-RAISED CONFCRS BLOCKAGE 
Cc (INPUT) 
C 
C NRAC —- NUMBER OF RACIAL PGSITICNS CN SURFACE 
C CINFLT) 
C NANG — NUMBER OF ANGULAR POSITICNS ON SURFACE 
€ (INPUT) 
c NPCINT — NUMBER GF ARBITRARY POINTS CN SURFACE 
c (NRAD¥NANG) 
C LJ — JET VELOCITY (INPUT) 
C LM — MAIN STREA® VELOCITY (INPUT) 
C C - ACTUAL DIAMETER OF JET CRIFICE (INPLTI) 
C R ~ MON-CIMENSICNAL RAOTAL PCSITION OF SURFACE 
C PINTS (GIMENSIONED AT LEAST NRAD) (INPUT) 
C THETA — ANGULAR POSITICN OF SURFACE PCINTS 
C (DIMENSIONED AT LEAST NANG) 
C xP —- X-PQOSITION CF SURFACE POINTS 
C (CIMENSTONED AT LEAST NECINT) 
C YP - Y-POSITICN CF SURFACE PCINIS 
C (DIMENSICNEC AT LEAST NPCINT) 
C x - HCRIZONTAL COCRDINATE OF JET 
C 2 - VERTICAL CCCRCINATE OF JET 
C ZINC — NQN-CIMENSICNAL INCREMENTAL STEP SIZE FCR Z 
C CELTAZ - 2-PCSITION CF EFFECTIVE/VIRTUAL SQURCE 
C CELCCK - EFFECTIVE BLCCKAGE CYLINCER CIAMETER 
€ RKSAV — CUMULATIVE VCRTICITY ALCNG THE JET 
c XI — NATURAL CCCROINATE TANGENT TO JET 
C iET£ — NATURAL CCCRCINATE NGRMAL TC JET 
C ETA - NATURAL COCRDINATE PERPENCICULAR TC 3I 
C ANO ZETA 
C ALPFA — ANGLE BETWEEN X-AXIS AND XI-AXIS 
C FAREA — INCREMENTAL PLATE AREA 
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C Teg EeR 


eee RE SCELTANT FORCE ON FLATE 
C XBAF X=-PGSITION CF RESULTANT FCFCE 
C YT NORMALIZEC INCUCED FCRCE 
PMEETCit RESL 48 (A-F,C-Z) 
PEN ae ibm CARN OO Vs TEETAC SO} 
CIMENSICN XP 9YF(45C0),2F(459) 
CATA INPUT 


AC(2,60) NRAD 
AC (2,60) NANG 


LT STCRAGE MATRICES 
T = NRAD*NANG 


UNO 


5C ,50 
(50,5 
(450) 


xO 
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C VELCCITY RATIO 
VR = LI/ZUM 

c fee riCN OF EFFECTIVE SCURCE 
CELTAZ = VR¥DSQRTI(VRI/15.0 

C CCRRECTICN FCR DELTAZ = 0.0 
IFQNTYPESLE1L) CELTAZ = 0.001 


WRITE(6,8C) DELTAZ 
WRITE(6,90) ZINC 
40 FORNAT(2F8.4) 
EC FORMAT(*0', 2X y*JET VELOCITY = %, F8.4s2X,"UNIFCRM', 
X*VELCCITY = '9F8.4,2X, "JET CIAMETER = ',F 8.4) 
EC FCRNAT(I4) 
iC FORMAT(F1C.5) 
gC FORMAT("'O%, 2X y* EFFECTIVE SOURCE IS %yF5.242X, 
X*CIAMETERS AGCVE THE PLATE') 
$C FORMAT (*°G%,2x,°Z STEP SIZE = 'yFS.2 
€ CEFINITION CF PI 





63 





Y) 
ult 
- 
<I 
ee 
tt 
© 
ow 
O 
© 
WO ——~ 
om om (1) OD) 
= HAN < ~~ 
qa 2921WNW m4 (NCL) 
- VIvpuIt <_1 cic 
~ WM WwuOoUY Lae 
y) 3 CJR 4 bt mt 
> — — r- ttm WWW) 
— Pc z. cCc7—or HHH 
one Lu oO I mm _) — om 
UO = WY CJ) jj ~w ie Ey ty 
O ~ tu mal wwe O aad add 
= WY) or Iw lb aaxx eo wee woe 
uw Oo w OO OW wWw O «xt f 20.0. 
~ > x» DD = ~N Fe i ee Sh MOK 
oO “sa — mm <{ WY etn x >< tee = 
oO 'o) ' > ud) — OO Ren OL wrt KC 
OO <— tts Cc) °° @ aa) uw) + + Sent Sane Soe? 
© iy _~ ox. + mw sy aa KN me e- ++ 
w on oOo Ww om tt) ~—— tt fren esANQD Ree eo qu (te 
Oo Lu = a ww ~ WwW [22 am ANN ade WOrYr WOrt &. wmnaa 
Oo th q Zz © ~~ tO) er we |— f-—- Om Mww S27 2tet 2S ee ASS s4aN@O 
> Ff or = wrt O Oo qIwnz NN NSN wT Tr dat OFO OFM ZZ. thet aI <d<f 
4 2 uw — ww Of ny 2 a 2 -t-Or Odt¢d AN NZZ FRR LLU Lub WY PHN MMW 
~ «A t -— oe. ~~ I CC <f mW) Ze ONIN me tet wee Ot tO CS OoOwU 
<— tt ww £ a ve «£ <-_ —_ i Le Mma rt ee COA MA FO OR ON CIV 
+ WwW kK He CD - <I ~* UW Y)?) Krew NU Fr NOG Set www www DOORN OO! % Ht Ht 
WY a F aI — fH FH OC WW wee er dd =~ ODO MW FENN FRET 2AZ NNW 2 MRAIN™N aN 
wm CC Cee uw © — wu Lom NZ { <f Oo oO + OO NOOO acct eee OOUO &§ NNw © AAS 
aq WU Pp OQ O <a jw LL " @ OmuUwWU c) ® ++ OW Nee QHD YVMM YOU O ww! | awl NINN 
# it cm Cc ke <I ONY a Oo Oo = OO we ADL OOO wv 13 we le ee 
oOo - WW 2&2 i ad “sa ™ Oa OO Nea) wae ee tt CQ nn a 
> “x aA O tK wook& 7) =~ Won ee n NNN #8 WHT ron wee CO wo >~ wut 
N YU \ 2 wae z OIw OO SHH OY NI mq woOodw W aN 
H~—- VY &- Vo W TO Ww FIO fe WY —_om = et Men 2862 SND etd WIM o&. W et (NIG 
tw ww <1 Ww Rts BS A Pe mrs Uhl tC he ana VWRew AQSNM NTS HN aad 
Ls 1 ot COC <a nus fF WW we ww wv) mN 22a UU 2tet VwVYW — UOrFmgmw << KrF 
_ “D =) 1daruy wo kh th Orm a tN YY ANAIN AN <aieteQ SU) Rete OUD Metre FR LU 
Oo 2. ©W WO ww VU OO VY WSR OU KR ONIN LONI AK ORR Re SES PUY WU SYS ROKK WONN™ 
cy —) ee < AL t 
t— “I us m Vv) <{ 
WO =z 
(a) 
oe) 
WYO YWoY OS 1 Vw ) mn | > | w 








SeerelwLe CONSTANTS 


C 


OR my OTD OT OTD ote HD oy 
Te tp TO 
MMMNUIMMMYW 


+oeHorer+ 


WILL tL OC 
OMNNNMNMWY) 


+e+ +++ + + 

(NO 
NAVE HANA Ie{IMNnoonw 
TUR HHEHEHHANNNNNNYNWUY) 
HH HAL FH FH we et ee ee 
ITV OCWLAINAt RR eee eee 
RRS Ede e{ Je CK CCS 
WWF OOF RF OF OOUOIOIU0OU 
md pet 4 et LILI LI NIM) Y) 
W WO >So OS NI RYINILULIOOOQOOW OW 


—Jeet ey ng hehe te eke te te te ee 
Az. 

Mr AO OO eat QDI Ore 
WLLIVIVIVIVIVIVIVIVIYVILICILICICIWOC)UI 


AaAMm 
bbe 
LLILLI 5€ d€ 


ag? Mee? es te 
wen 


mtNCD AT 
a i 
WUWM 
UUVUY 


ZETAI*CCAL 


Ll 
L 


ZETA2*CCR2 


2 


ZETAB*CCA3 


ZETAB*CONS 


V4 


XICNE#COAL 


= 
>_> 


wi 
h2 
h3 
4 


XITWC*CON2 


—= 
—= 


IMI EC CN: 


ET PEYECN 4 


MENTAL VORTICITY TERM 


— 
S 


INCE 


OO 


C((DEXP(Z2/BI-DEXP(Z1L/BIIS¢C1L CHLEXP((Z24+Z1)/E))) 
(1 .0/8.40)4# VR*¥2 * DATAN(FF) 
HK + HKSEZV 


FF = 
+ K 


FK 


HK 


EK SAV 
INCREMENTAL INTERFERENCE VELOCITY COMPONENTS 


C 





C.0 
_ pe eee re tuasecnmnenne a) SSHNAGSUNSCE SEE) 


wh 

f 

= 

‘ 

mt 

tN -> 

ww <I ws 

wl) #t 

Ow 

or 

Oo 

(NI 

we tf e tf 
maw 

aa } a... 

od ad * 
xs 

UO WU 


© 
ma 
= 
O 
= - 
a 
=> WwW 
~~ © 
Ww 
oOo = 
a < - 
a QO Oo 
© 
++ W + © 
Y oO 
== 7 em OO 
oe her ee ee ° 
ee OW D OQ 
me OF fw ® 
ww tl ww UW 
a>—lUCc rlCUrMCi 
= «7 ° 
wnO oO F- 
WwW Y) 
-—— l wu) 
“> © \ ca 
v7 Ww re 
1 + LL. ts) 
aay wa L$ iw 
aI rr eS 
wv) 
iJ 
= 
Ww 


65 





EG iC 100 
{10 Sen TiNue 
meeoriLCUL4T ION OF BLOCKAGE VELOCITY TERS 

TFONTYPE.~EQ.C) GO TC 120 

IFCNTYPE.EQ.2) GO TC 120 

RELCCK = C-96sCSQRT (VR) 

C8LCCK = 


“Ceo D 
PrP OOK 


ea th Wt 


b—4 p=4 


Geet CULATION CF PRESSURE CCEFFICIENT 
120 CPhI 55) = ~4eC¥H(UCT SI) + UCT, J)*F2 + VII, I)442) 
CFANGE CESTCNATICN OF VARTAELE 
ZP(C) = CF(I,J) 
28 SNTINYE 


C 
C 
C 
i N 
1 N 
C CALCLLATION OF INCREMENTAL FORCE 
Q 
¢ 


ay ada’ “1 TT) P> 2) 
Pr KMOMZOPAeO) 
BD *SV*ACDADM 
—a— mm Mmneoamnn it i 


- 


RAGE PRESSURE COEFFICIENT 


CPAVG = (CP(I,J) + CPIK,JII/2.0 
memecGLETICN OF EFFECTIVE FCRCE 
FORCE = PAREAACPAVG 
Cie sto C 


D 

EQ.1) 

EQeNANG) GO_TC 150 
= 2.0*FCRCE 


66 





arte t Ri )/ 2e0)=CEOS(IBETAC J) } 


CALCLULATICQN 


XM 


CF MCMENT ABOUT Y-AXIS 


WU 


FORCE*XM 


YMC 
SUMMATICN CF TOTAL FORCE 


C 


TOTFCR + FCRCE 


ICTFCR 
SUMMATION CF TOTAL MOMENT 





Cc 


TYMCM + YMC 


—_ 
_ 


TYMCM 


BeecCUtevT IGN CF CENTER CF ACTION 


> 


TYMCM/ TOTFOR 
CALCLLATICN CF NCRMALIZEC INDUCED FORCE: Y/T 


XE AF 


e 


1.0/VF 
2eO*TCTFCR*VRIAV*¥*2/PT 


CATA CLTPUT 


VRINY 
YT 





C 


~ 
- 
— 
- 
~ 
x 
VW" x 
rt 
Ys — re 
a. MN ™ 
UO e e 
- oO ff 
° LL. LL. 
as —~ & ve 
Ov — x >< 
tad Ty o~t 
BS “> & o~ 
ae weer (NY - 
= = @ 
- -~O WL 
»~< am iL © 
a — 
~ ~ x tL) 
AN = YO 
e ——~ & a 
sO TN © 
iL ~ @ 5 
— a om © 
Qa - mi oF 
~ ~* = 
© " “>< <I 
e wer (NJ = 
oS] <j au a 
Oo UN IN 
mt UL ~ @ Y) @ 
O# JIT OQ~ 0d UO 
ZwenOr AOL CC UL. 
dm = (Y=w © - © 
ar x1 OM nm 
oc >< ~ mm aS 
att OQ tee ON Ni ft 
wI©D @« QO o ~ 
WLoe toue - x< 
t— +1 O NI QO 
fo ow mm ™ Se WII fj 


tt} Owes Ow NDw et 

oO wk [2 Owk 2k 
Nutu<l ew mettcrtrei<ity) 
INP SOONER > RP eee 
Oma «Mm wre ever 
UVZaAUnUYCONIUIOMWVA 
Waste Oo suVUOowe 
x x< 


OOVO 
O MOT) 
INE 


yc 


© 
oO 
a) 


Oo O20 VY 
Tr wo ~™ 
Mme INI tu 


67 








Or 


eye 


i 


3. 


Pee. Ror RE CES 


English, J. W., "The Design and Performance of Lateral 
Thrust Units for Ships, Hydrodynamic Considerations," 
leds. R-E.N.A., VOL. 105, 1963. 


English, J. W., and Steele, B. N., "The Performance of 
Lateral Thrust Units for Ships as Affected by Forward 
Speed and Proximity of a Wall," Ship Division, N.P.L., 
mem. NO. 28, 1962. 


Stuntz, G. R., and Taylor, R. J., "Some Asvects of Bow- 
Meruscer Design,” Trans. S.N.A.M.E., Vol. 72, 1964. 


Chislett, M. S. and Bjorheden, 0O., “Influence of Ship 
Speed on the Effectiveness of a Lateral-Thrust Unit," 
Hydro-Og Aerodynamisk Laboratorium, Rep. No. Hy-8, 
April 1966. 


Principles of Naval Architecture, S.N.A.M.E. 


Bradbury, L. J. S. and Wood, M. N., “The Static Pressure 
Distribution Around A Circular Jet Exhausting Normally 
from a Plane Wall into an Airstream, C.P. No. 822, 

mmot., A.R.C., 1965. 


Margason, R. J., "The Path of a Jet Directed at Large 
Angles to a Subsonic Free Stream," NASA TN D-4919, 1968. 


Gordier, R. L., "Studies on Fluid Jets Discharging 
Normally into Moving Liquid," St. Anthony Falls Hyd. 
Lab., Tech. Paver, No. 28, Series B, Aug 1959. 


Sucec, J. and Bowley, W. W., "Prediction of the 
Trajectory of a Turbulent Jet Injected into a Cross- 
flowing Stream," Journal of Fluids Engineering, Trans. 
A.S.M.E., Series I 98. Dec 1976. 


Abramovich, G. N., The Theory of Turbulent Jets, MIT Press, 
1963. 


Jordinson, R., "Flow in a Jet Directed Normal to the Wind," 
nemeno, 3074, Brit. A.R.C., 1958. 


Keffer, J. F. and Baines, W. D., "The Round Turbulent Jet 
a Crosswind,” Journal of Fluid Mechanics, Vol. 15, 1963. 


Volger, R. D., "Surface Pressure Distribution Induced on 
a Flat Plate by a Cold Air Jet Issuing Perpendicularly 
from the Plate and Normal to a Low-Speed Free-Stream 
Flow," NASA TN D-1629, 1963. 


68 








14. 


Dis 


Ol 


17. 


Nore 


JSR 


20. 


Zee 


Ae 


aS. 


Fearn, R. L. and Weston, R. P., "Induced Pressure 
Distribution of a Jet in a Crossflow," NASA TN D=-7916, 
LoD. 


McMahon, H. M. and Mosher, K.D., “Experimental Investiga- 
tion of Pressures Induced on a Flat Plate by a Jet 
Issuing into a Subsonic Crosswind," NASA SP-218. 


Kamotani, Y. and Greber, I., "Experiments on a Turbulent 
Jet in a Crossflow," A.I.A.A. Journal, Dec 1974. 


Wooler, P. T., Burghard, G. H. and Gallagher, J. T., 
"Pressure Distribution on a Rectangular Wing with a 
Jet Exhausting Normally into an Airstream," Journal 
G@emAtreratt, Vol. 4, 1967. 


Wooler, P. T., "On the Flow Past a Circular Jet Exhausting 
at Right Angles from a Flat Plate or Wing,” Journal of 
the Royal Aeronautical Soeiety, Vol. 71, March 1967. 


Wu, J. C. and Wright, M. A., "A Blockage-Sink Represen- 
tation of Jet Interference Effects for Noncircular Jet 
Orifices," NASA SP=-218. 


Adler, D. and Baron, A., "Prediction of a Three Dimensional 
Circular Turbulent Jet in Crossflow," Research Report, 
Israel Institute of Technology, Haifa, Israel, Seotember 
1978. 


Schmitt, H., "Anderung einer Parallel stromung entlang 
einer ebenen Platte durch einen quer gerichteten 
Freistrahl," Zeitschrift fur Flugwissenschaften und 
Weltraum forschung, Sept/Oct 1979. 


Durand, W.F., Aerodynamic Theory, Vol. II, SECT l, 
Dover 1963. 


Robertson, J. M., Hydrodynamics in Theory and Application, 
Prentice Hall Inc. 1965. 


69 





PieeeAG DESTRIBULION LIST 
No. Copies 


Defense Technical Information Center 2 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0142 2 
Naval Postgraduate School 
Monterey, California 93940 


Department Chairman, Code 69 Mx i 
Department of Mechanical Engineering 

Naval Postgraduate School 

Monterey, California 93940 


Professor R.H. Nunn, Code 69 Nn 5 
Department of Mechanical Engineering 

Naval Postgraduate School 

Monterey, California 93940 


LCDR Bradford B. Waterman 1 


468 Oak Street 
East Bridgewater, Massachusetts 02333 


70 








bait 90180 


W229968 Waterman 
Cel Analysis of jet~ 
crossflow interactions 


with application to 
Ship bow thrusters, 








aie 
a 





i 


$ 
01 
NOX LIBRARY 


| 


ii 


3 


| 


9 
2 
U 





, 
_ 
; 





‘3 
OT ag 
; 7 
os 
BLEW 


re 





‘ o a, 


«a 





